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Abstract 

In a homogenous and isotropic cosmology, we introduce general exact solutions for 
some modified gravity models. In particular, we introduce exact solutions for power- 
law f(R) gravity and Brans-Dicke theory in Einstein and Jordan conformal frames. In 
the Brans-Dicke case, the solutions are presented for both single and double exponential 
potentials in Einstein frame which correspond to power-law potentials in Jordan frame. 
Our analysis for extracting general exact solutions can also be generalized to those scalar- 
tensor theories in which the scalar field has an exponential coupling to Ricci scalar. 
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1 Introduction 

General Relativity is a powerful tool to explain theoretically many observational facts about 
the Universe. Despite all the successes, there are also some unresolved problems such as in¬ 
flation, the cosmological constant problem and the problems associated with the dark sector, 
i.e., dark matter and dark energy. These problems have motivated people to seek for some 
modifications of the theory. Among many possibilities, much attention has been paid in recent 
years to two classes of modified gravity theories; scalar-tensor theories [1] in which gravity is 
described not only by a metric tensor but also by a scalar field and f(R) gravity theories [2] 
in which R in the gravitational action is replaced by the function f(R). 

Comparing with General Relativity, these theories consider new degrees of freedom which open 
new possibilities for addressing the aforementioned problems. In scalar-tensor gravity, one in¬ 
troduces a new dynamical scalar field and in f(R) gravity one considers forth order derivatives 
of the metric in the held equations. These modifications make the theories admit a larger 
variety of solutions than Einstein equations in General Relativity at the cost of increased com¬ 
plexity. This increased complexity together with nonlinearity make hireling exact solutions for 
the held equations be much more difficult. 
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Despite the complexities, some cosmological exact solutions have been found in both classes. 
In scalar-tensor gravity, exact solutions are usually given when the scalar field has no potential 
function [1], There are also exact solutions for some functional forms of the potential such as 
power-law potential [3] and exponential potential when coupling of the scalar field is minimal 
[4], In power-law f(R ) gravity, exact solutions have been reported in homogenous and isotropic 
cosmology in [5] [6]. In the present work we are looking for exact solutions in the two modified 
gravity models with a different approach. We start with Einstein frame representation of a 
general scalar-tensor gravity in which there is an interaction between the scalar field and mat¬ 
ter systems. This representation is the same for all f(R) gravities and all parameterizations 
of scalar-tensor theories. Different gravitational theories are characterized by their potential, 
their coupling functions and the coupling strengths (A(<p) and /3((p) in the following). We 
introduce general exact solutions of the gravitational models for which the coupling strength 
takes constant values. 

The plan of the paper is the following : In section 2, we first consider the action of a general 
scalar-tensor theory in Einstein frame in which the matter system is taken to be a perfect 
fluid with a barotropic equation of state. We then write the corresponding field equations in 
a homogenous and isotropic cosmology describing by the Friedman-Robertson-walker (FRW) 
metric. We assume spacetimes with zero spatial geometries. The matter system and the scalar 
field are not separately conserved due to the interaction and there is an energy flow between 
the two components. Direction of the energy flow depends on the coupling strength and the 
equation of state of the fluid. For an equation of state corresponding to radiation there is 
no interaction and hence no energy flow. In section 3, we present the solutions of the field 
equations by direct integration. We first solve the continuity equation for the matter system. 
The solution indicates that evolution of the matter energy density is modified due to the inter¬ 
action. We can consider this modification as a modification of the exponent of the scale factor 
by an arbitrary function e. Here e measures the energy exchange between the matter and the 
scalar field. This leads to a relation between the scalar field and the scale factor. A relevant 
role is played by this relation to reduce the second order differential equation of the scalar field 
to a first order one for obtaining the potential function. We will show that exact solutions 
can be obtained when e and the coupling strength take constant values and for exponential 
forms of the potential. We argue that this procedure can be applied to power-law f(R) gravity 
and some scalar-tensor gravities such as BD theory. In section 4, we draw our conclusions. 


2 Field Equations 

We consider the action functional 

Sef = U dV=0 {MpR - <r - 2V (93)} + s m (A 2 (<p) g ^, VO ( 1 ) 

where M ~ 2 = 87 tG, G is the gravitational constant and S m is the action of some matter field 
-0. The function A(<p) is a coupling function that characterizes coupling of the scalar field p 
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with the matter sector. The action (1) may be taken as Einstein frame representation of a 
generalized scalar-tensor gravity or an f(R) gravity model . 

Variation of this action with respect to the metric tensor leads to 


g„„ = + T”) 


- flU J 


where 


T - -g^pV.p - V (<p)g l 


fits 


rjirn = ~2 SSmig^tp) 


( 2 ) 

( 3 ) 

( 4 ) 


are stress-tensors of the scalar field and the matter system. Variation with respect to the scalar 
held p, gives 


□<£> 


dV{v) ^ /%) , 

dp 


where 


m = m, 




dlnA(p) 


dp 


( 5 ) 


( 6 ) 


and T m = g^ v T™. The two stress-tensors T™ and T^ v are not separately conserved. Instead 
they satisfy the following equations 


= -V M 7X = ^^-V v p T n 


- fllS 


- /J,V 


M„ 


( 7 ) 


We apply the held equations in a spatially hat homogeneous and isotropic cosmology described 
by FRW spacetime 

ds 2 = — d,t 2 + a 2 (t)(dx 2 + dy 2 + dz 2 ) (8) 

where a(t) is the scale factor. To do this, we take T™ as the stress-tensor of a perfect fluid 
with energy density p m and pressure p m . In this case, the gravitational equations (2) takes the 
form 


3— — M p ( p m + p^) 


( 9 ) 


a a 


2 I y ~ (Pm T Pip) 


a a 


( 10 ) 


where p v ^ p 2 + V (p) , p v = -j p 2 — V(p) and overdot indicates differentiation with respect 

to the cosmic time t. Combining these equations, gives 


a a" , , 9r l 

f —^ 


+ ^^2 ““ [~(Pm — Pm) + V ((/9)] 


From (5) and (7), we obtain 


.+3 

a dp 


P(<P) 

M n 


(pm 3 Pr 


( 11 ) 


( 12 ) 
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where 


Pm + 3— (u m + 1 )p m — Q 
a 

(13) 

Pcj> + 3—(cp^ + 1 )p lf = —Q 
a 

(14) 

n ■( q \ 

Q = M VKPm 3 p m ) 

(15) 


is the interaction term, uj m = p m jp m and u v = p^/p v - The direction of energy transfer depends 
on the sign of Q. For Q > 0, the energy transfer is from the scalar held (dark energy!) to the 
matter system and for Q < 0 the reverse is true. 


3 Solutions 


There are three independent held equations among (9)-(14) for finding a(t), (p(t) and p m (f)h 
The equation (13) can be easily solved 


Pm 


p0m & 


—3(cj m +l) 


(1 —3u >m) 


f fidtp 


(16) 


in which pom is an integration constant. This solution indicates that the evolution of energy 
density is modified due to interaction of ip with matter. There will be no loss of generality if 
we write (16) as 

Pm = POm. a _3(u,m+1 ) +e (17) 

with e being defined by 

(1-3 u m )f/3d(p 

€ = - M^Ta - (18) 

Even though e and the coupling strength [3 in (18) are generally evolving functions, we will 
restrict ourselves to the case that they can be regarded as constant parameters. In this case, 
(18) reduces to 

<p — aM p In a (19) 

with a being a constant defined by the relation e = /5cr(1 — 3co m ). The expression (17) is similar 
to the rule presented by some authors for characterizing decaying law of vacuum energy into 
dark matter [8]. It states that when e > 0, matter is created and energy is constantly injecting 
into the matter so that the latter will dilute more slowly compared to its standard evolution 
p m oc a _3 ( u,m+ i). Similarly, when e < 0 the reverse is true, namely that matter is annihilated 
and the direction of energy transfer is outside of the matter system so that the rate of dilution 
is faster than the standard one. 

We will show that (19) actually satisfies (12) for some potentials. It is possible to find those 

^The action (1) and the subsequent field equations are similar to those of some interacting models in which 
there is an interaction between dark energy and (dark) matter [7]. In those models, the coupling function A(ip) 
is usually taken as a pre-assumed function that is specified, for instance, by phenomenological arguments. 

iHere we assume that the perfect fluid that describes the matter system has a constant equation of state 
parameter u) m . 
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potentials by solving a first order differential equation. There are two important cases that 
the function /3((p) takes a constant configuration; f(R) gravity and Brans-Dicke theory. We 
consider the two cases separately in the following : 


3.1 f(R ) Gravity 

The action for an f(R ) gravity theory in Jordan frame is given by 

S JF = l -M 2 p I d 4 xV^f(R) + S m (g^ (20) 

where g^ is the metric in Jordan frame. We consider a conformal transformation 

9nv = A ~ 2 { { f) 9nu ( 21 ) 

with A~ 2 (ip) = ^ = f'(R). This together with 

■f ^jlnA(p) ( 22 ) 


and /? 



transforms (20) into the action (1) with a potential 


[9] [10] 


V(<p{R)) 


M 2 R 
2 


m 

f'\R) 


(23) 


Instead of solving the field equations for a given f(R) function (or, equivalently, a given 
potential function), we are looking for those models that accept the solution (19). In this way, 
one can reduce the second order differential equation of (p to a first order one for finding the 
functional form of V(tp). To do this, we first put (19) into (12) which leads to 


a a 

~ + 2— = 
a cr 


<7 A/2 


P , , . 1 dV(<p) 

{fim 3Pm) 17 M p d<p 


(24) 


Comparing the latter with (11) gives then a consistency relation 

dV(<p) a 

+ 17 ^ (V 7 ) — a Pm 


d^p 


M n 


(25) 


where a = ^— 1) + T?-(3a; m — 1). This together with (17) and (19) gives 


dV(<p) 

dip 


a 

+ K 


V(<p) = ap 0m e l ~ 3{uJm+1)+£],p/(TMp 


This is a first order differential equation which leads to the following solution 


(26) 


V(ip) = + 5p 0me \rU^+i)+,\vl7M, 


(27) 
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where 


( 28 ) 


^ _ § \Wm ~ 1 ) + / 5(3 Ul m ~ 1 ) 

[— 3(uj m + 1 ) + fda{l — 3u m )]^ + a 

and C\ is an integration constant. Thus the expression (19) is an exact solution of the held 
equations for the potential (27). For C\ = 0, (27) is reduced to a single exponential potential 


V{<p) = Sp 0m e [ - 3{oJm+1)+e]ip/aMp 


(29) 


Putting this potential into the Friedman equation (9), gives 


— = u a 
a 2 


3(cJ m + l)+€ 


(30) 


where 

POm ^ + 1 

V = - 

3 — cr 2 /2 

For a 2 6 , this integrates to 

a ( f ) = {7(vV t + c 2 )Y h 


(31) 

(32) 


with 7 = l[3(cj m + 1) — e]. For e < 0, 7 remains positive for a perfect fluid which satisfies 
weak energy condition uj m + 1 > 0. However when e > 0, there is a chance for 7 < 0 and a 
contracting scale factor. This case happens when the energy how from p to matter is so strong 
that make the Universe collapse. 

To answer the question of how this solution is attributed to f(R) function, we should use a 
back mapping of our results to the Jordan frame. One may use (19), (21) and (22) to write 


a(t) = A(<p)a(t) = e M pa(t) = a /?<T+ 1 (f) = { 7 (y/u t + U 2 )} ( - /3cr+1 ^' /7 (33) 

di = A(lp) dt = dt = a / 3 o "(f) dt = { 7 (y/u t + C , 2 )} /3ct ^ 7 dt 
=>t = {7 (y/u t + c 2 )}^ +1 /\/u(/3a + 7 ) 

Combining the latter two relations gives d(t), scale factor in the Jordan frame 
scale factor has a power-law form both in Einstein and Jordan frames : 

1 _ /3<t+1 

if co m = 1/3 =>■ a(t) oc 1 2 , a(t) oc t f )a + 2 

2 _ 2(/3<t+1) 

if uj m = 0 =>■ a(t) oc 13-/30- , a(t) oc t t/^+s) (35) 

These solutions indicate that even hough passing from one conformal frame to the other does 
not change the power-law behavior, it changes the exponents. However, in the case that 
/da « 1, the exponents of the scale factor follow the standard evolution in the radiation and 
the dust phases in both conformal frames. 

§The case C\ A 0 corresponds to a double exponential potential which will be considered later. 


(34) 

. For C 2 = 0 , 
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The function f(R) depends crucially on the functional form of the potential V(ip). For a single 
exponential potential of the form (29), the expression (23) gives a differential equation 


V(R) = Sp 0m f'J 


Mp.R f 


(36) 


This has a simple power-law solution f(R) = fo R n with n being 

l)((3cr + 1 ) 


( 


n = 


[Wm + 1) + /5cr (cu m + 1/3) 


and f 0 is given by 


25 


POm. 
M2 

p 


K nf 0 




(37) 


(38) 


Thus (32)-(35) are exact solutions for power-law f(R) gravity model. It should also be noted 
that n is given in terms of two free parameters u m and a (or e) so that our solutions consider 
a much wider variety of solutions comparing with those obtained in [6]. 


3.2 Scalar-Tensor Gravity 

The general action of a scalar-tensor gravity is given by [11] 

SjF = J {F(cf>)R - Z{(t>)r v y^ V u (j) - 2 U(<j>)} + S m (9pu, tl> m ) (39) 

where F((j>), Z((j)) and U((f>) are some functions^. This action is reduced to the action (1) by 
the conformal transformation ( 21 ) with A{ip) = F - 1 / 2 ( 0 ) and 

(*£)* = 2 M 2 \-( dlnm ) 2 + 2^-1 (40) 

Vlp 4 df> )+ 2 F(0) J 1 j 

V(<p) = M 2 f/( 0 )F~ 2 ( 0 ) (41) 

The coupling function A(ip) depends on the functions F((/)), Z{4>) and U (0) through the relation 
(40). For some particular choices of these functions, (3 takes a constant configuration and then, 
as a result of ( 6 ), A(ip) takes an exponential form. These cases define a class of scalar-tensor 
theories for which the solution (19) and all the subsequent results obtained in the subsection 3.1 
are valid. In the following, we will restrict ourselves to this class of scalar-tensor theories. One 
important theory in this class is given by the BD parameterization in which F(«p) = 16nG(f), 
Z(4>) = IQfGubd/4> and = 87 tGW( 0 ), and then 

S JF = f d^x^g^R - - W{<!>)) + S m {g^ 0) (42) 

with ix bd and W{4>) being BD parameter and the potential in Jordan frame, respectively. This 
action is reduced to ( 1 ) by [ 1 ] [ 12 ] 

A(<p) = e ^ BD ' p/Mp (43) 

^One can always redefine the scalar field to reduce F((f>) and Z((j)) to only one unknown function. 
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( 44 ) 


(p(<j>)/M p = Jujbd + 3/2 ln(-^-) 

0o 

V((p) = W(<f>(<p)) e 8/W/Mp (45) 

where 0 O ~ G and /3bd = + 3/2. The exact solutions of this BD model with 

potential (29) are given by (32) just by replacing (3 with /3bd- Note that when cj B d 0, then 
f3 bd ~and Einstein frame representations of BD model and /(/?,) gravity are the same. 
In Jordan frame, on the other hand, the solutions are (33)-(35) with 0 being replaced by [3bd 
and the corresponding potential is of power-law form 

= Sp 0m (^f + ^ ] (46) 

00 

The solutions are plotted in hg.l for both f(R) and BD gravity models. The figure indicates 
that the distinction between the two theories becomes more apparent when time grows. 

There is a large lower bound imposed on ujbd by solar system experiments [13] which implies 



Figure 1: The plot of scale factor for dust in (a) Einstein and (b) Jordan frames. The solid 
and dashed lines correspond to a = 0.5 while dotted and dashed-dotted lines correspond to 
cr = 1 in f(R) and BD gravity, respectively. 


that (3bd « 1- If cr remains of order of unity, then /Jg^cr « 1 and the exponents of the scale 
factor in the power-law solutions (35) follow the standard evolution in radiation and matter 
phases both in Einstein and Jordan conformal frames. 

The condition for existing a late-time accelerating phase is given by q — —da/a 2 < 0 with q 
being the deceleration parameter. For the solution (32) in the dust phase, this condition gives 
cr > l/f3 (in the BD case a > 1/(3bd)- Since (3bd < (3 then accelerating solutions exist for 
larger values of the parameter a in the BD theory with respect to f(R) gravity. 

In potential (27), one can take di / 0 which leads to a double exponential potential. The 
integration constant C\ can be determined by noting the fact that when ip takes a constant 
configuration in the action (1), then V(tp) acts as a constant cosmological term. One may 























consider the condition 


(47) 


V(<p — 0) — C\ + dpam - A A = 

which results in C\ = p\ — 8p 0m . The relation (27) takes then the form 

V(<p)/p 0m = (pa/P om - 8)e-^ M r + <fe[ -Z(“m+l)+eM*M p ( 48 ) 


For this potential function, the Friedmann equation (9) becomes 


— = u a 
a 1 




+ v a ^ 


(49) 


where 


PmO Pa/P mO 8 

Ml 3 — cr 2 /2 


(50) 


For a 2 6 , exact solutions can be found by a direct integration for different values of uj m . We 
first consider some popular cases. 

1 ) 1 • 

In this case, (28), (31) and (50) give 8 = —1, u = 0 and v = 3 ( 3 -^71) w h ere we have set 

Pa/P om = 7/3 according to recent observations [14], The solution of (49) is then 


a 


(t) ~ t + C 3 )] 


2/cr 2 


(51) 


where C 3 is an integration constant. The deceleration parameter is negative if cr 2 < 2 and the 
solution leads to a power-law inflation for a sufficiently small cr. 

2) u m = 1/3 : 

In this case, 8 = 12 ff 3a 2 , u = ^f 3(4 / fT2) and v = We have simple solutions 

in some special cases 

a(«) = [(( + C 4 ) 2 -| 1/2 if u 2 — 2 (52) 

a(t) = [4^2/3 t + 2 C 5 } 1/2 if a 2 = 7/2 (53) 

The case a 2 = 7/2 and C' 5 = 0 corresponds to the standard Friedmann model a(t ) oc f 1//2 in 
the radiation phase. The solution for any a is 


cr 2 (ua a2 


+ va 4 ) 


{2a 


3 +< t ‘ 


U 


+ va 


2 —< 7 2 . 


'1 + 


ua 


a 2 — 4 


- 2 F 1 I 


cr 


> 1 + 


cr 


ua 


2-4 


2’ 2(cr 2 — 4) ’ 2(cr 2 — 4) ’ 


-]} 


— t+Ce (54) 

where 2Fl[a,b,c,x] is the hypergeometric function 2 F\ [a, 6 ; c; x]. Note that this expression 
gives t = t(a). 

3) 0 . 

In this case, 8 = /J'^// 3 /l[^ a and € B d = Pbd& and the solution of (49) is 


o 2 (ua a (P BD+ ^ + va 3 ) 


{2a 2+Cy2 UcfiBDcr-l + va 2 -v 2 ^ 


ua~ 3+ ^ BDa+cj2 r l 

1 +-2F1 


a " 


v 


2 ’ 2(—3 + (3 BD cr + a 2 ) ’ 
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<7- 


ua 


-3+/3 BD (T+ct 2 


2(—3 + 0bd& + & 2 ) 


-]} — t+C 7 


(55) 

This solution is plotted in fig.2 for some values of the parameter a and Cj = 0. The general 



0 2 4 6 8 10 12 14 


a 


Figure 2: The Plot of t — t(a ) in the dust region (the solution (55) for C 7 = 0). The solid 
lines correspond to BD theory for /3bd = —0.08 (or ujbd = 4000) and a = 0.5, 1 from top 
to bottom. The dashed lines indicate f(R) gravity (/3 = —1/\/6) for the same values of the 
parameter a. 


solution of (49) for any cr and oj m is given by 

-1 

[e B D ~ 3 (U m + 1 )}[ua a2+e BD + va 3(o, m +l)] 


{2a CT +( 2 + 3w ™) yj + +Va ~<T 2 +2 



-2F1[-, 

u 2 


£bd — 3 (uj m + 1) 

2 [cr 2 + 6bD — 3 (u!m + 1)] 


£bd — 3 (a j m + 1) 
2[a 2 + 6bd — 3 (u!m + 1)] 


yd~ 


u 


]} — t + c 8 


(56) 


where 6 bd = @bdg{}- ~ 3u; m ). It should be remarked that the solutions (51)-(56) can also be 
regarded as the solutions of Einstein frame representation of an f(R) gravity with potential 
function (48) just by replacing /3bd and 6 bd with (3 and e, respectively. The f(R) function 
which corresponds to this double exponential potential can be formally obtained by using 
(23) although it is not easy to solve the resulting differential equation for general a and 
U! m . However, it is shown that [15] the Einstein frame representation of the models^ f(R) = 
R — /i 2 t n+1 ) / R n has a potential function which takes a double exponential form when n —> oo. 
This potential corresponds to (27) in some regions of the parameters space. 

The fig.2 plots the solution (55) for the two gravity models. It indicates that the overall 
behavior of the solutions is the same for given values of the parameter cr. It is important 
to note that there is no distinction between the two gravity models in the vacuum and the 


II These models have been investigated in [16] to explain the late-time acceleration of the Universe. 
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radiation phases since the solutions (51)-(54) do not depend on the parameter (3. 

We may consider different scalar-tensor gravity models for which (6) takes a constant value 
and then the coupling function A[p) becomes exponential. Besides BD theory, in all these 
cases (19) satisfies the held equations for some potential functions of the forms (29) and (48). 
For instance, consider F(0) = F 0 e v ^ and Z(0) = Z 0 e v ^ for which (40) gives p = \f2C,M v <\) 
with C = y/fh 2 + Then (6) results in the coupling function A(p) = F 0 x l 2 e ^IM p in 
which /3 C = is a constant. This parameterization has been already studied as extended 

quintessence model where the scalar held playing the role of the dark energy is exponentially 
coupled to the Ricci scalar [17]. For the potentials (29) and (48), this scalar-tensor theory has 
general exact solutions of the forms (32)-(35) and (51)-(56) in which (/?, e) and (Pbd, £bd) 
should be replaced by ({3 C , e c ). 

4 Conclusion 

We have introduced homogenous and isotropic exact cosmological solutions in some modified 
gravitational theories. We assume that the Universe has a hat spatial geometry and is filled 
with a perfect fluid with a barotropic equation of state. In a general scalar-tensor theory in 
Einstein frame, there is a coupling between the scalar held and the matter part through the 
coupling function A(p). We have found general exact solutions for gravitational theories in 
which the coupling strength fi(p) takes constant values. There are two groups of gravitational 
theories with a constant coupling strength, f(R) gravity and a class of scalar-tensor gravity 
including BD theory. Exact solutions of these gravitational theories are found for single and 
double exponential potentials. 

For single exponential potentials, the exact solutions are hrst presented in the Einstein frame 
for f(R) gravity and BD theory and then transformed back the solutions to the Jordan frame. 
It is shown that the solutions in Jordan frame belong to a power-law f(R) gravity and BD 
theory with a power-law potential. Analysis of the solutions indicates that there are power-law 
solutions for the scale factor in Einstein and Jordan conformal frames for both gravitational 
theories. When f3cr « 1 (/3bd& « 1) they reduce to the standard evolution in radiation and 
dust phases. 

For double exponential potentials, general exact solutions are also presented for both BD and 
f(R) theories in Einstein frame. It is interesting that these solutions do not distinguish be¬ 
tween the two theories in vacuum and radiation phases. The analysis can be generalized to 
scalar-tensor theories in which the scalar field has an exponential coupling to Ricci scalar. In 
this case, functional form of the solutions are exactly the same as the solutions in the other 
two gravitational theories for single and double exponential potentials and can be obtained by 
usual replacement of the parameters with (f3 c , e c ). 
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